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The diurnal variation of hydrogen, nitrogen, and chlorine 
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Abstract. In situ measurements of hydrogen, nitrogen, and 
chlorine radicals obtained through sunrise and sunset in die lower 
stratosphere during SPADE are compared to results from a 
photochemical model constrained by observed concentrations of 
radical precursors and environmental conditions. Models allowing 
for heterogeneous hydrolysis of N 2 0 5 on sulfate aerosols agree 
with measured concentrations of NO, N0 2 , and CIO throughout 
the day, but fail to account for high concentrations of OH and H0 2 
observed near sunrise and sunset. The morning burst of [OHJ and 
[H0 2 ] coincides with the rise of [NO] from photolysis of N0 2 , 
suggesting a new source of HO, that photolyzes in the near UV 
(350 to 400 nm) spectral region. A model that allows for the 
heterogeneous production of HN0 2 results in an excellent 
simulation of the diurnal variations of [OH] and [H0 2 ]. 

Introduction 

Data collected during the Stratospheric Photochemistry, Aero- 
sols, and Dynamics Expedition (SPADE) provide a unique oppor- 
tunity to study the diurnal variation of hydrogen, nitrogen, and 
chlorine radicals in the lower stratosphere. Flight tracks were 
designed, using a high resolution meteorological forecast model, 
to sample air with nearly uniform concentrations of precursors of 
free radicals through sunrise and sunset on successive days. In 


harvard University 

2 NOAA Aeronomy Laboratory 

3 Jet Propulsion Laboratory 

4 NASA Ames Research Center 

5 NASA Langley Research Center 

6 University of California 

’University of Denver 

8 National Center for Atmospheric Research 

9 NOAA Climate Monitoring and Diagnostics Laboratory 

l0 Goddard Space Flight Center 

Copyright 1994 by the American Geophysical Union. 

Copyright 1994 by the American Geophysical Union. 

Paper number 94GL02782 
0094-8534/94/94GL-02782S03 .00 


this paper we compute the 24 hour variation of [OH], [H0 2 ], 
[NO], [N0 2 ], and [CIO] using a photochemical model constrained 
with measured values of pressure, temperature, aerosol surface 
area, concentrations of radical precursors (0 3 , H 2 0, CH 4 , NO , 
and Cl, inferred from CC1 3 F and C0 2 ), total column ozone, and 
planetary albedo. Comparison of calculated and observed concen- 
trations of radicals at different times throughout a daily cycle pro- 
vides a rigorous test of our current understanding of photochemi- 
cal processes that regulate the abundance of radicals. 

Theory and Observation 

We focus on flights of the ER-2 from Moffett Field, CA 
(37.4°N) on 11 May 1993 (930511) and 12 May 1993 (930512) at 
"“19 km altitude (-65 mbar). Concentrations of N 2 0 and 0 3 were 
nearly uniform from before sunrise to mid-morning on 930511 
and from mid-afternoon through sunset on 930512 (Figure 1 ). Air 
parcels with values of [N 2 0] between 240 and 260 parts per bil- 
lion (ppb) and pressure between 64.5 and 69.5 mbar were selected 
for analysis in this study. 

Inputs to the photochemical model are given in Table L Con- 
centrations correspond to the mean of measured values for 
selected air parcels, except for [NO,] and [CH 4 ] which were 
estimated using correlations with [N 2 0] from other flights during 
SPADE. Concentrations of inorganic chlorine (Cl,) were com- 
puted from in situ measurements of halogenated source gases and 
[C0 2 ] [Woodbridge et al., 1994]. Total column ozone and ultra- 
violet albedo were specified from observations by the Total Ozone 
Mapping Spectrometer (TOMS). 

Reaction rates and absorption cross sections from DcMore et 
al. [1992] were used, except for extrapolation to longer 
wavelengths of cross sections for H 2 0 2 , HN 0 2l and HN 0 4 . A 
reaction probability of 0.1 was used for the heterogeneous hydro- 
lysis of N 2 0 5 [DeMore et al., 1992], The heterogeneous hydro- 
lysis of CINO 3 was included, but has a negligible effect on model 
results at the temperature (214.5 K) of these observations. Photo- 
lysis rates were calculated using a radiative transfer model that 
accounts for Rayleigh and aerosol scattering [Prather, 1981] using 
vertical profiles for [ 0 3 ] from a climatology scaled to match the 
total column measured by TOMS, and vertical profiles for tem- 
perature and aerosol extinction from the National Meteorological 
Center and the Stratospheric Aerosol and Gas Experiment II, 
respectively. 

We assumed that each species reached a balance between pro- 
duction and loss, over 24 hours, for the temperature, pressure, and 
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Figure 1. Measurements of pressure* [N 2 0], and [0 3 ] obtained 
during ER-2 flights 930511 and 930512 plotted as a function of 
solar zenith angle ("am" and "pm" distinguish observations 
obtained before and after local noon, respectively). Dark points 
indicate air parcels, with pressure between 64.5 and 69.5 mbar and 
[N 2 0] between 240 and 260 ppbv, selected for analysis. 

latitude at which an air parcel was sampled. Since winds were 
weak during May, the flow was primarily zonal, and photochemi- 
cal time constants arc short, we expect on average relatively small 
deviations from steady state. Further details arc provided by 
Salawitch et al. [1994]. 

Model results and observations arc shown in Figure 2. Calcu- 
lations neglecting heterogeneous chemistry underestimate [OH], 
[H0 2 ], and [CIO] and overestimate [NO] and [N0 2 ] by large 
amounts. For the lowest solar zenith angles (SZA), a model 
allowing for the heterogeneous hydrolysis of N 2 0 5 on sulfate 
aerosols ( JPL-llet , Figure 2) simulates the concentration of all 
radicals to within the uncertainty of measurement (2o measure- 
ment uncertainties, including random and systematic effects, are 
estimated to be 30% for [OH], 40% for [H0 2 ], 20% for [NO], 
25% for [N0 2 ], and 25% for [CIO]). However, this model fails to 
account for the early morning rise and late afternoon decay 
observed for HO* ([OH]+[H0 2 ]) [Wennbcrg et al., 1994]. 

The lime dependence of the model ‘.measurement ratios for 
[NO], [CIO] and [OH] is illustrated in Figure 3. To evaluate the 
significance of these comparisons, we assume: 

[XL™ = S [X]mi;as +Z + P (1) 

where [X] A ™ and [XImiias represent the actual and measured con- 
centrations of species X, and S, Z and P represent the systematic 


Table 1. Model inputs 


Latitude 

37.9°N 

o 3 

1 .4 ppmv 

Temperature 

214.5 K 

h 2 o 

4.4 ppmv 

Pressure 

66.9 mbar 

ch 4 

1 .5 ppmv 

Column 0 3 

353 DU 

NO, 

5.0 ppbv 

UV Reflectivity 
Surface Area 

0.34 

6.0x10 8 cm' 1 

Cl, 

1 .2 ppbv 


error, zero artifact, and precision of the measurement, respec- 
tively. Nighttime observations of radicals that disappear in the 
dark can be used to estimate Z. The short term variability of the 
observations, given the nearly uniform concentration of precur- 
sors, can be used to estimate P. The ratios in Figure 3 are res- 
tricted to times when measured concentrations exceed Z+P by a 
factor of 3. In this case the ratios should differ from unity by a 
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Figure 2. Measurements (dots) of [OH], [H0 2 ], [NO], [N0 2 ] 
(crosses and dots represent data from the JPL and NOAA instru- 
ments, respectively), and [CIO] during ER-2 flights of 930511 and 
930512, for selected air parcels illustrated in Figure 1, plotted as a 
function of solar zenith angle. Model simulations using rates and 
cross sections from DcMorc et al. [1992]: neglecting all hetero- 
geneous processes ( JPL-Gas , blue dotted line); including hetero- 
geneous hydrolysis of N 2 0 5 and C1N0 3 ( JPL-llet , blue solid 
line); and same as JPL-llet except allowing for updated absorption 
cross sections for HN0 3 and quantum yields of 0( l D) and the 
heterogeneous production of HN0 2 (Model B\ red line). Measure- 
ments of [CIO] and [NO] (obtained every 16 and 1 seconds, 
respectively) have been smoothed using a 1 minute median filter. 
[NO], [N0 2 ], and [CIO] have been divided by [NO,] and [Cl,] to 
account for slight variations in the abundance of radical precur- 
sors. Values of [NO]/[NO y ] calculated using JPL-Gas (-0.18 at 
mid-day) have been omitted to emphasize comparisons for the 
heterogeneous models. 
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Figure 3. The ratio A(X) s ([X] M0D ,; L -[X| MKAS )/[Xj M0DI , Lt where 
[X]model and [X]mhas represent modeled and measured concen- 
trations of species X, for the JPL-llet (circles) and Model B 
(crosses) simulations. Ratios have been averaged into SZA bins. 
The ratio is displayed only when [X] MKAS > 3 x(Z+/ 5 ) (see text). 


constant amount, 5, if the model describes properly Urc response 
of radical concentrations to changes in solar illumination. 

Differences between calculated and observed [NO] and [CIO] 
during early morning are close to the uncertainty of measurement, 
and show no trend with SZA (Figure 3), demonstrating that the 
variation of photolysis rates for N0 2 and C1N0 3 with solar illumi- 
nation is calculated reasonably well. It is unclear if the trend with 
SZA between calculated and observed [NO] during late afternoon 
is real, or what process is responsible. However, the trend with 
SZA of the model-measurement ratio for [OH] (Figure 3) exceeds 
the measurement uncertainty by large amounts, illustrating the 
existence of a source of HO* during early morning and late after- 
noon not accounted for by the JPL-llet model. 

The largest sources of HO* at low SZA are reactions of 0( ] D) 
with H 2 0 and CH 4 , and photolysis of HN0 3 . Peak production of 
0( J D) from photolysis of 0 3 and peak photolysis of HN0 3 occur 
at wavelengths near 305 nm; at high SZA, solar radiation for this 
spectra] region does not penetrate to ER-2 altitudes. The abrupt- 
ness of the morning rise of [HO*] and its coincidence with the rise 
of [NO] suggests that the missing source of HO* is due to rapid 
photolysis of a gas at wavelengths longer than 320 nm, the spec- 
tral region of peak photolysis of N0 2 . Significant production of 
0( 1 D) at these wavelengths can be ruled out based on the analysis 
of Michel sen et al. [1 994] . 

Wennberg et al. [1994] suggested that photolysis of HN0 2 
(which is rapid and peaks at wavelengths near 370 nm) or a 
molecule with similar photolytic behavior is responsible for the 
missing source of HO*. To match the morning rise of [HO*] 
requires [HN0 2 ] = 10 parts per trillion (pptv). However, [HN0 2 ] 
calculated for sunrise using the JPL-llet model (the primary sinks 
and sources of HN0 2 are photolysis and the three body reaction of 


OH and NO, respectively) is too low, 0.04 pptv, to provide an 
appreciable source of HO*. 

Heterogeneous processes have been suggested to explain pro- 
duction of HN0 2 in the troposphere [Notholt et al., 1992] and 
environmental chambers [Pitts et ah, 1985]. Most of these studies 
have focused on reactions involving NO* and H z O. Production of 
HN0 2 in the stratosphere by a heterogeneous reaction involving 
N0 2 could explain the rise of [HO*] at sunrise, following photo- 
lysis of HN0 2 formed throughout the night. However, Lhis pro- 
cess would have to be slow to limit production of HN0 2 during 
the night to ~10 pptv, and therefore would be unable to produce 
enough HN0 2 during the day to account for the discrepancy 
between modeled and measured [HO*] at late afternoon. Produc- 
tion of HN0 2 by a rapid reaction of a non-radical species with an 
abundance comparable to [HO*] is required to explain both the 
morning and afternoon discrepancies. Two possible candidates, 
H 2 0 2 and CH 2 0, have production rates too slow to provide appre- 
ciable amounts of HO* during late afternoon. 

Wennberg et al. [1994] proposed heterogeneous decomposi- 
tion of HN0 4 , 

HN0 4 (aerosol)— > HN0 2 + 0 2 , (2) 

as a possible source of HN0 2 sufficient to account for the source 
of HO* at high SZA. In aqueous solutions HN0 4 undergoes uni- 
molecular dissociation producing HN0 2 and 0 2 [Logagcr and 
Sehested, 1993; Kcnlcy et al., 1981]. The slow rate constant 
measured in these studies, for a temperature of 300 K, suggests 
this mechanism is loo slow to explain the atmospheric 
observations unless the Henry’s Law coefficient for HN0 4 is 
extremely large. HN0 4 has been observed to stick to ice surfaces 
with high efficiency, but without becoming chemically active [Z. 
Li, R. R. Fricdl, and S. P. Sander, in preparation, 1994]. However, 
Zhu et al. [1994] suggested heterogeneous decomposition of 
HN0 4 was the source of HN() 2 in their laboratory chamber. We 
are unaware of any studies of the chemistry of HN0 4 on sulfuric 
acid. 

We examine the consequence of the heterogeneous production 
of HN0 2 using Model B , identical to the JPL-llet model, except 
allowing for: heterogeneous decomposition of HN0 4 on sulfate 
aerosols to form HN0 2 with a reaction probability of 0.2; quantum 
yields of 0( ] D) from photolysis of 0 3 -25% larger than recom- 
mended values [Michelsen et al., 1994]; and measurements of the 
temperature-dependent absorption cross section for HN0 3 [Burk- 
holder et al., 1993] that result in a 20% reduction in the photolysis 
rate of HN0 3 at 220 K. Different photolysis rates for 0 3 and 
HN0 3 used in Model B result in a decrease in calculated [NO*], an 
increase in [CIO], and negligible change in calculated [HO*] since 
variations induced by each change cancel. This model results in 
an excellent simulation of the latitudinal variation of radicals 
observed during SPADE [Salawitch et al., 1994]. 

Agreement between theory and observation for [OH] and 
[H0 2 ] is improved dramatically using Model B (Figures 2 and 3). 
Calculated concentrations of HN0 2 at night (12 pptv) are 
sufficient to produce a large source of HO* at sunrise. The supply 
of HO* from photolysis of HN0 2 is almost exhausted prior to the 
onset of significant production from 0(*D), resulting in a local 
maximum in [HO*] during early morning. The rale of decay of 
[HO*] during late afternoon is reduced because HO* that is con- 
verted to HN0 4 at midday is recycled by reaction (2) and subse- 
quent photolysis of HN0 2 . 

Measurements of [HN0 4 ] and [HN0 2 ] in the lower strato- 
sphere would greatly aid our understanding of processes govern- 
ing the diurnal variation of [HO*]. Profiles of [HN0 4 ] have been 
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obtained in the lower stratosphere from a weak Q-branch absorp- 
tion feature near 803 cm -1 recorded by the Atmospheric Trace 
Molecule Spectroscopy Experiment during May, 1985 [Russell et 
al., 1988]. However, measured [HN0 4 ] has large uncertainties 
due to limited laboratory data for the temperature dependence of 
the absorption coefficient of HN0 4 and is consistent with values 
calculated by models that both include or ignore reaction (2), 
which results in nearly a factor of 2 reduction in calculated 
[HN0 4 ] for 20 km, at sunrise. We are unaware of any stratos- 
pheric measurements of [HN0 2 ]. 

Conclusions 

Models that include heterogeneous hydrolysis of N 2 O s on sul- 
fate aerosols provide an excellent description of the variation with 
solar illumination of observed concentrations of nitrogen oxide 
and chlorine radicals in the lower stratosphere. A model that 
allows also for the heterogeneous production of HN0 2 reproduces 
key features of [HO,] observed during SPADE: the burst at sun- 
rise coincident with the appearance of elevated [NO], the local 
maximum during early morning, and lire slow afternoon decay. 
Observations of the diurnal variation of radicals during SPADE 
shown here suggest that rates for chemical loss of 0 3 by cycles 
involving HO, are probably more rapid than previously believed 
[Wcnnberg et al., 1994]. 
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